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ABSTRACT 

Using deep HST ACS imaging and VLT F0RS2 spectra, we determined the velocity dispersions, ef- 
fective radii and surface brightnesses for four early-type galaxies in the z = 1.237 cluster RDCS 1252.9- 
2927. All four galaxies are massive, > IQ^^Mq. These four galaxies, combined with three from 
RDCS 0848-1-4453 at z = 1.276, establish the Fundamental Plane of massive early-type cluster galax- 
ies dX z = 1.25. The offset of the Fundamental Plane shows that the luminosity evolution in rest-frame 
B is AlnM/is — (— 0.98±0.06)z for galaxies with M > 10^^-^ Mq. To reproduce the observed mass- 
to-light ratio (M/L) evolution, we determine the characteristic age of the stars in these M > 



galaxies to be 3.0_o;3 Gyrs, i.e. z* = 3.4_q;4. Including selection effects caused by morphological bias 

(the "progenitor bias"), we estimate an age of 2.ltQ 2 Gyrs, or = 2.3to;2 for the elliptical galaxy 
population. Massive cluster early-type galaxies appear to have a large fraction of stars that formed 
early in the history of the universe. However, there is a large scatter in the derived M/L values, which 
is confirmed by the spread in the galaxies' colors. Two lower mass galaxies in our z = 1.25 sample 
have much lower M/L values, implying significant star-formation close to the epoch of observation. 
Thus, even in the centers of massive clusters, there appears to have been significant star formation in 
some massive, M ~ 10^^ Mq, galaxies at z ~ 1.5. 

Subject headings: galaxies: clusters: general — galaxies: elliptical and lenticular, cD, — galaxies: 
evolution — galaxies: fundamental parameters — galaxies: photometry — galaxies: 
clusters: RDCS 1252.9-2927 
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1. INTRODUCTION 



The Fundamental Plane (FP) allows a direct measure 
of the mass and the mass-to-light ratio, M/L, of early- 
type galaxies. The FP combines three variables, the ef- 
fective radius (ve), the average surface brightness within 
the effective radius (/g), and the velocity dispersion (a 



These data are combined i nto the relation a^-^^ oc Tel^ 
for the rest-frame B band l|,Torgensen et alJll996h . With 
such quantities, we can measure M/L oc u'^/ijele) and 
how it depends on mass, oc TeCr^. Massive galaxies out to 
z ~ 1 appear to ev c )lve as A In M/L b ::::: z for both clus- 
ters llKelsoii et all l2nnnbl: Ivan Dokkum fc Franxl WM. 
IWuvts et a l."2nn4') and in some field samples, though 
there is a larger scatter f orthe latter (van Dokkum ct aL 
2^01"; Gcbhardt et al. 2003; van Dokkum & Ellis ■200'^ 
va,n de Ven et al...200,'^ .van der Wei et aL 2004b^. This 
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slow rate of evolution implies an early epoch of forma- 
tion, z/ ~ 3, for the stars in early- type galaxies assum- 
ing passively evolving simple stellar populations. How- 
ever, at z ~ 1.25, only ~ 50% of stellar mass we ob- 
serve today have been formed (e.q.\ 'Mad au et a'nil998l: 
Stcidcl ct al. 1999; Rudnick ct al. 2003). This impHes 
that the majority of stars in cluster early-type galax- 
ies formed long before the average star in the universe. 
Observations determining the luminosity-weighted age of 
galaxies close to z = 1.25 will test this, and there are only 
three gal axies with FP measurements t o date at these 
redshifts ijvan Dokkum fc Stanforll2003|) . 

We observed four luminous early-type galaxies in the 
z — 1.237 rich, massive and X-ray luminous c luster of 
galaxies RDCS 1252.9-2927 ijRosati et all 120041 ) using a 
combination of the Very Large Telescope (VLT) Focal 
Reducer/low dispersion Spectrograph 2 (F0RS2) and 
the Advanced Camera for Surveys (ACS) on the Hubble 
Space Telescope (HST). RDCS 1252.9-2927 is the most 
massive cluster found to date at z > 1, thus it contains 
a number of luminous and, likely, massive galaxies. We 
use the FP to constrain the M/L evolution and set the 
mass scale for four gala xies in RDCS 1252.9-29 27. These 
results, combined with Ivan Dokkum fc Stanfo rd ( 2003), 
measure the ages of the stellar populations in early-type 
galaxies at z = 1.25. We assume a Vim = 0.3, f^A = 0.7 
and Ho — 70 km s~^ Mpc~^. All observed magnitudes 
are in the AB system. However, for comparison with 
previous work, we convert observed magnitudes into rest- 
frame Johnson B using the Vega zeropoint. 

2. RDCS 1252.9-2927 FUNDAMENTAL PLANE DATA 

Four galaxies in RDCS 1252.9-2927 we re selected from 
among the nine known cluster members ijDemarc n i2nul 



2 



iLidman et al.ll200^ with zgso < 22.5 mag that fit into a 
single multi-sht mask. These are the first, second, third 
and fifth brightest cluster members. An image of each is 
shown in Figure ^ along with the spectra, in descending 
order of brightness. Below we discuss the measurement 
of Te and le from the ACS data, the a from the F0RS2 
spectra and the conversion to the rest-frame Johnson B. 
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Fig. 1. — Mosaic showing, from left to right, the observed spectra 
near 4000 A in the rest-frame of the cluster, images and images of 
the residuals. Each image is 2'.'0 on a side with 0'.'05 pixels from 
the 2g5o ACS data, corresponding to roughly rest-frame Johnson 
B. The fit results are listed in Tabled The absolute value of flux 
in the residuals is <10% of the flux in the original data. For #4419 
and #4420, the residuals have an "S" shape which is interpreted 
as a sign of interaction fBlakeslee et al. 20Q3)- The residuals for 
#6106 show an over subtraction in the center, it is likely that the 
r^^'^ law is too steep; in this case the best fitting profile is oc r^^'"^. 



2.1. HST ACS Imaging 

As described in iBlakeslee et all ()200 3D , the Advanced 
Camera for Surveys (ACS) imaged RDCS 1252.9-2927 
with four overlapping pointings. Each pointing has three 
orbits in the F775W filter, or 1775, and five orbits in the 
F850LP filter, or ^§50. A de Vaucouleurs, or r^/^ model, 
was fit to each of the individual z^^n images using the 
method of lvan Dokkum fc Franxl l)1996(l . Each galaxy in 
each image had a unique point spread function generated 
using the TinyTIM v6.2 package (Krist 1995). The two 
galaxies at the middle of the cluster, referred to as # 4419 
and # 4420 in Tabled were fit simultaneously. Tabled 
contains the average of the best fitting parameters for 
each image with the galaxies listed in order of decreasing 
zgso flux. We plot in Figure ^ a mean image, corrected 
for the ACS distortion, of all the zsso images for each 
galaxy along with the average residuals from the fits. 

The product rg/g is used to measure the evolution in 
M/L. Because of the strong anti-correlation between the 
error for /ig and the error for (Joreensen ct al. 1993), 
the uncertainties on this product for all four galaxies is 
small at ~5%. 

2.2. VLT F0RS2 Spectra 

The four galaxies in Table ^ were observed using 
F0RS2, on the VLT, through slit masks with the 600z 
grism in conjunction with the OG590 order separation 
filter. The observations were done in service mode with 
a series of exposures, dithered over four positions, for a 
total integration time of 24 hours. The resulting signal- 
to-noise (S/N) ratios at 4100 A rest-frame are listed in 



Table Details concerning the dat a reduction are de- 
scribed in Ivan der Wei et all Il2004ar . 

The high spectral resolution, 80 km s~^ per pixel, re- 
sulted in accurate internal velocity dispersions for the 
four cluster member s (see Table [B . The s p ectra were 
fit, by the met hod oflvan Dqkkuni fc Franxl l)1996il . with 
stellar spectra ijValdes et al.ll2004D with a wide range in 
spectral type and metallicity. For more details concern- 
ing the usage of the t emplates and the derivatio n of ve- 
locity dispersions, see Ivan der Wei et all ll2004ar . 

The velocity dispersions were aperture corrected to a 
1.7 kpc ci rcular aperture a t the d istance of Coma as de- 
scribed in I.Torgensen et iiD 1)19961) . The listed errors in- 
clude a statistical error derived from the value of the 
fit and a systematic error estimated to be at most 5% for 
the spectra with the lowest S/N ratio. 

2.3. Rest-frame Magnitudes 

In order to compare with other FP results in the liter- 
ature, the observed zgso magnitudes must be converted 
into Brest , the equivalent of observing the galaxies with 
a rest-frame Johnson B filter (Bcssell 1990) in the Vega 
system. The ^§50 filter is centered at 4058A in the rest- 
frame of the galaxies in RDCS 1252.9-2927. This filter is 
close to the central wavelength of Johnson B at 4350A, 
but even the modest wavelength difference means that 
the conversion between zgso sd z — 1.237 and the John- 
son B will depend on the color of the galaxy. To compute 
this conversion, we redshif ted the Sbc and E templates 
from IColeman et alJ l|198nD to z = 1.237 and calculated 
the Brest magnitude as a function the observed zgso and 
i775~J color, yielding Brest = zsso-0. 45(1775- J)-|-l. 68. 
This approach is sl igh tly different than that used i n 
IKelson et a1.l ^mO) orlvau Dokkum Sta,nfordl p003|). 
but yields results that differ in the mean b y < 0.02 mag- 
nitud es, with an error of only < 5% (see iHolden et alJ 
120041) . The J, along with K^, photometry comes from 
th e VLT ISAAC and N TT SOFI observations discussed 
in ILidman et al.l l)2004|) . We will also use this data to 
examine the rest-frame optical B — I colors below. All 
colors were measured within an aperture of two effec- 
tive radii. The ACS imaging was smoothed to match the 
seeing in the ISAAC data to measure this color. 

The statistical errors on the FP are dominated by the 
error on the velocity dispersions, which are around 10% 
including an estimate of the systematic error. Because 
this error dominates the error budget, we will take the 
error on a to be the FP error for the rest of the paper. 

3. THE FUNDAMENTAL PLANE AT Z = 1.25 

The most straightforward way to measure the evolu- 
tion in M/Lb is to compute the offsets for the seven 
gala xies at z = 1.25 from the FP of the Coma clus- 
ter l|Jorgensen et al.lll996j) . In Figure |21 we plot our 
data along, with the FP for Coma, using the rest- 
frame B. We show the average offset for the five 
most massive galaxies out of the total of seven in our 
sample. We will discuss below why we remove those 
two lower mass galaxies. Measuring the offset from 
the FP ensures that we measuring AM /Lb for galax- 
ies at the same part of the FP or, roughly the same 
mass. The offset AM / L and rate of M/Lb evolution 
is readily apparent in Figure |3| We find AM /Lb = 



The Cluster Fundamental Plane at z = 1.25 
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Fig. 2. — A projection of the FP for our sample and 
Coma. The results for RDCS 1252.9-2927 {z=1.237) are plotted 
as solid red dots, the green crosses are from RDCS 0848+4453 
^an Dokkum & Stanford 2003, z = 1.276) , and the blue open 
stars are from Coma I Jorgens en et al.lll996l z = 0.023). The y- 
axis is xj-p = 0.83 log re + 0.69 log /e, one of the projections of the 
FP. The solid line is the FP for Coma while the dotted line has 
the same slope but is shifted AM /Lb = —0.98 for the five massive 
galaxies at z = 1.25. 



— 1.23 ± 0.08 for five z = 1.25 early-type galaxies, three 
from RDCS 1252.9-2927 and two from RDCS 0848H-4453 
with masses M > 10^^'^ Af©. This corresponds to an 
evolution in AIuM/Lb oc (-0.98 ± 0.06) z, a small 
deviation from the AlnM/£ R oc (-1.06 ± 0.09) z of 
Ivan Dokkum fc Stanford (2003) and the AIuM/Lb oc 
-1.08z of lWuvt s et al. ( 2004). 

There is a large scatter seen in FigurcElin AlnM/Ls, 
(T{\n M/Lb) = 0.32 for the seven z = 1.25, early-type 
galaxies. This scatter is twice the size of the scatter in 
Coma or MS 1358+62, regardless of whether the scatter 
is computed for all galaxies, or only the seven most lumi- 
nous galaxies in either MS 1358+62, or the Coma clus- 
ter sample. A large part of this scatter comes from the 
two lower mass galaxies in the sample. The five galaxies 
with M > 1O"-^M0 show cr(lnAf/Ls) = 0.22, which is 
not statistically different from the Coma or MS 1358+62 
value. There is an obvious selection effect towards low 
M/L galaxies in a luminosity-selected sample. This may 
both increase the scatter and bias the mean change in the 
M/L ratio, hence we remove the two low mass galaxies 
from our sample. 

The M/Lb values for the M > 10"-^Mq galaxies at 
z = 1.25 correlate with the rest-frame B — I colors, 
as seen in Figure 0| Both the colors and the M/Lb 
track a rapidly declining star-formation rate model from 
l^ruzual & Chariot (2003). As this relatively simple stel- 
lar population reproduces most of the observations, the 
observed scatter in M/Lb is then likely the result of a 
spread in the luminosity-weighted ages. 

At lower redshifts, the M/L is a function of the total 
mass (Jorgensen ct al. 1996). Such a trend is observed at 
z = 1.25 (see Figure El where lower mass galaxies have 
lower M/Lb and, therefore, bluer colors (see Figure 01) 
The slope of the mass-M/L relation appears to steepen 
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Fig. 3. — Change in M/Lb for early-types as a function of 
redshift. We use the same symbols as Figure 1^ with the addition 
of light blue open squares for MS 1358+62 I K elson et al.l l2000bt 
z = 0.328), grey open "+" are from MS 2053-04 {z = 0.583) and 
the fi lled orange tr i angle s are from MS 1054-03 (z = 0.832), both 
from IWuvts et ahl 120041^ . The resulting evolution, with respect 
to the Coma FP, is shown as a line with the form AlnAf/L^ <^ 
(-0.98 + 0.06)2. 



Coma 




Fig. 4. — M/Lb as a function of the rest-frame Vega B — I 
color for galaxies with M > 10^^'^Mq. The average Coma val- 
ues are represented as a blue star, while the three galaxies from 
RDCS 1252.9-2927 are solid circles and the two galaxyes from 
RDCS 0848+4453 are green crosses. The line shows the trajec- 
tory of a solar metallicity model with an exponentially declin- 
ing, T = 200 Myrs, star-formation rate Bruzual & Chariot (200^ 
model normalized to the average Coma cluster observations. The 
highest M/Lb is is #4419, the brightest cluster galaxy. The col- 
ors of #4419 are significantly bluer than predicted for its observed 
M/Lb, ruling out the offset being a result of dust or metallicity 
effects. The lowest mass galaxy in RDCS 0848+4453 is excluded 
(see text). 



at higher redshifts. This is expected for a population 
where the spread in the M/L comes from the spread in 
age. As a stellar population becomes younger, the M/L 
changes more quickly, so the spread in M/L will grow 
as observations probe closer to the epoch of formation. 
However, this trend will be exaggerated by the selection 
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Fig. 5. — M/ Lb as a function of mass, using tiie same symbols as 
Figure 121 but not including the results for MS 2053-04. The lower 
redshift trend of higher M/Lg at higher masses appears to be 
preserved at z = 1.25, even when ignoring the lowest mass galaxy 
in RDCS 0848+4453. The dotted lines represent Lb = lO" Lq 
and 10^° Lq. Our 2850 = 22.5 selection limit corresponds to is ~ 
IO^^'^Lq for a galaxy with the colors of a Coma elliptical. 



of galaxies at a fixed luminosity as the sample selection 
will prefer lower M / L galaxies. Lower mass galaxies with 
larger values for M/Lb will not appear in this sample 
because of our magnitude limit, effectively Lb — IO^^Lq. 

4. IMPLICATIONS FOR THE EVOLUTION OF 
EARLY- TYPE CLUSTER GALAXIES 

The rate of the M/ L evolution can be used to constrain 
the luminosity- weighted age for early- type galaxies. The 
most straight-forward estimate is to assume all galaxies 
formed at one epoch and find the age of the galaxies that 
will produce the observed change in M/L from z — 0.023 
to z = 1.25 using population synthesis models. For our 
sample of five M > galaxies, the mean age is 

T* = S.O^o'a Gyrs before the time of the observations, or 



a formation redshift of ; 



3.4l°;4 using the same mod- 



els as fvan Dokkum fc StanfordI l)2003D . namelv IWorthevI 
(^^). This agrees with the results from the high mass 
sample of[Wuvts et al. ( 2004), who found = 2.95 
for galaxies regardless of morphology. 



+0.81 
0.46 



When computing the age for the early-type galaxy pop- 
ulation, there is an overestimate of the age of a popula- 
tion caused by young galaxies not being counted as part 
of the early-type population, even if those young galaxies 
will evolve into early-types after the epoch of observation. 
This "progenitor bias" depresses the rate of observed 
evolu tion in M/L by up to 20% ijvan Dokkum fc Franxl 
1200 1|) . Using this same assumption, namely that the true 
evolution is AIuM/Lb oc (—1.18 ± 0.06) 2;, we instead 
Gyrs before z = 1.25, or a forma- 



find 



= 2.1 



+0.2 



■^-0.2 

tion redshift of = 2.3lo2. iBlakeslee etall l|2003D and 
Ekiman ct al. ( 2004) both find a mean age of > 2.6 Gyrs 
using the color s of th e galaxies in RDCS 1252.9-2927. 
IBlakeslee et alJ l)2003ft re moved the "p r ogenit or bias" 
with simulations, whereas iLidman et alJ l|2004D uses all 
galaxies, regardless of morphology to similar effect. 

The above results imply that the stars that formed 
these massive galaxies were created at redshifts of ~ 
2 — 3, at which time less th an 1/3 of today's observed 
stellar mass was formed fe.o.: lBelill2005j) . However, there 
is a large spread in the M/Lb for all of the early- type 
galaxies at z ~ 1.25, larger than at lower redshifts. Us- 
ing colors confirms the spread, which can be interpreted 
as an underlying spread in the age of the populations. 
Thus, though ~ 2.5 for the most massive galaxies, 
some early-typegalaxies show much lower M/L values 
and, corresponding younger ages. In fact, the lowest 
M/L g alaxy in Figure 13) was tentatively classified by 
Ivan Dokk um fc Stanford ( 20i3) as having a recent star- 
burst based on the spectrum. Such younger appearing 
galaxies have lower masses than the high M/L galax- 
ies in our sample, but are still massive galaxies with 
log M/M0 ~ 11. The implication of all these results is 
that a significant fraction of the stars in the most massive 
galaxies appear to have formed very early in the history 
of the universe, before the majority of stars present to- 
day. The massive cluster galaxies appear to follow the 
same low-redshift trend of the higher mass systems hav- 
ing higher M/Lb- However, the larger spread at z = 1.25 
in the M/Lb indicates that we have identified some mas- 
sive, M ~ IO^^Mq, galaxies whose last burst of star for- 
mation occurred in the relatively recent past, z ~ 1.5. 

ACS was developed under NASA contract NAS5- 
32865, and this research was supported by NASA grant 
NAG5-7697. BH would hke to thank Daniel Kelson for 
useful discussions on the Fundamental Plane. The au- 
thors would also like to thank the referee for many useful 
suggestions. 
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TABLE 1 

DE VAUCOULEURS MODEL PARAMETERS AND VELOCITY DISPERSION 



Galaxy 












log Mass 


log M/Lb 




(") 


(mag per □") 




(km s-l) 


mag 


logMg 




4419 


2.806 


24.899 


24 


302 ± 24 


2.09 ± 0.02 


12.40 ± 0.08 


0.81 ± 0.06 


6106 


0.487 


21.573 


57 


294 ± 10 


2.07 ± 0.02 


11.61 ± 0.08 


0.22 ± 0.03 


4420 


1.016 


23.279 


29 


323 ± 21 


2.11 ± 0.02 


12.01 ± 0.06 


0.66 ± 0.06 


9077 


1.008 


23.529 


24 


130 ± 14 


1.90 ± 0.02 


11.22 ± 0.08 


0.01 ± 0.09 



^AU magnitudes arc AB. 

''Spectra have a resolution of 3.7 A(FWHM). 
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